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Summary. We review the main proposals of particle physics for the composition of 
the cold dark matter in the universe. Strong axion contribution to cold dark matter 
is not favored if the Peccei-Quinn field emerges with non-zero value at the end of 
inflation and the inflationary scale is superheavy since, under these circumstances, 
it leads to unacceptably large isocurvature perturbations. The lightest neutralino 
is the most popular candidate constituent of cold dark matter. Its relic abundance 
in the constrained minimal supersymmetric standard model can be reduced to ac- 
ceptable values by pole annihilation of neutralinos or neutralino-stau coannihilation. 
Axinos can also contribute to cold dark matter provided that the reheat temperature 
is adequately low. Gravitinos can constitute the cold dark matter only in limited 
regions of the parameter space. We present a supersymmetric grand unified model 
leading to violation of Yukawa unification and, thus, allowing an acceptable 6-quark 
mass within the constrained minimal supersymmetric standard model with /i > 0. 
The model possesses a wide range of parameters consistent with the data on the 
cold dark matter abundance as well as other phenomenological constraints. Also, it 
leads to a new version of shifted hybrid inflation. 



1 Introduction 

The recent measurements of the Wilkinson microwave anisotropy probe 
(WMAP) satellite on the cosmic microwave background radiation (CMBR) 
have shown that the matter abundance in the universe is Q m h 2 = 0.135^q'qqq, 
where Qi — Pi/p c with pi being the energy density of the i-th species and 
p c the critical energy density of the universe and h is the present value of 
the Hubble parameter in units of 100 km sec -1 Mpc -1 . The baryon abun- 
dance is also found by these measurements to be Hth 2 = 0.0224 ± 0.0009. 
Consequently, the cold dark matter (CDM) abundance in the universe is 
f2 C BMh 2 = 0.1126±gSi The 95% confidence level (c.l.) range of this quan- 
tity is then, roughly, 0.095 < J?cdm/i 2 ~ 0.13. Taking h ~ 0.72, which is its 
best-fit value from the Hubble space telescope (HST) [2], and assuming that 
the total energy density of the universe is very close to its critical energy den- 



2 George Lazarides 



sity (i.e. f2 tot ~ 1), as implied by WMAP, we conclude that about 22% of the 
energy density of the present universe consists of CDM. 

The question then is, what the nature, origin, and composition of this 
important component of our universe is. Particle physics provides us with a 
number of candidate particles out of which CDM can be made. These particles 
appear naturally in various particle physics frameworks for reasons completely 
independent from CDM considerations and are, certainly, not invented for the 
sole purpose of explaining the presence of CDM in the universe. 

The basic properties that such a candidate particle must satisfy are the 
following: (i) it must be stable or very long-lived, which can be achieved by 
an appropriate symmetry, (ii) it should be electrically and color neutral, as 
implied by astrophysical constraints on exotic relics (like anomalous nuclei), 
but can be interacting weakly, and (iii) it has to be non-relativistic, which is 
usually guaranteed by assuming that it is adequately massive, although even 
very light particles such as axions can be non-relativistic for different reasons. 
So, what we need as constituent of CDM is a weakly interacting massive 
particle. There are several possibilities, but we will concentrate here on the 
major particle physics candidates which are the axion, the lightest neutralino, 
the axino, and the gravitino (for other candidates, see e.g. Ref. [3 ). Note that 
the last three particles exist only in supersymmetric (SUSY) theories. 

In Sec. |U we examine the possibility that the axions are constituents of 
CDM. Sec.|3is devoted to outlining the salient features of the minimal super- 
symmetric standard model (MSSM), which will be used as a basic frame for 
discussing SUSY CDM. In Sec. 0] we summarize the calculation of the relic 
abundance of the lightest neutralino, which is normally the lightest supersym- 
metric particle (LSP), and investigate the circumstances under which it can 
account for the CDM in the universe. In Sees. [S] and H3 we discuss, respec- 
tively, axinos and gravitinos as constituents of CDM. In Sec. we present 
a SUSY grand unified theory (GUT) model which solves the bottom-quark 
mass problem by naturally and modestly violating the exact unification of 
the third generation Yukawa couplings. We study the parameter space of the 
model which is allowed by neutralino dark matter considerations as well as 
some other phenomenological constraints. Finally, in Sec. |H1 we summarize 
our conclusions. 



2 Axions 

The most natural solution to the strong CP problem (i.e. the apparent absence 
of CP violation in strong interactions implied by the experimental bound on 
the electric dipole moment of the neutron) is the one provided by a Peccei- 
Quinn (PQ) symmetry This is a global U(l) symmetry, U(1)pq, which 
carries QCD anomalies and is spontaneously broken at a scale f a , the so-called 
axion decay constant or simply PQ scale. Astrophysical (Sj and cosmological 
constraints imply that 10 9 GeV < / a ~ 10 12 GeV. The upper bound originates 
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|SJ El from the requirement that the relic energy density of axions does not 
overdose the universe. It should be noted, however, that this upper bound 
can be considerably relaxed if the axions are diluted (7J |SJ |S] by entropy 
generation after their production at the QCD phase transition (for more recent 
applications of this dilution mechanism, see e.g. Ref. [TU]). 

The axion is a pseudo Nambu-Goldstone boson corresponding to the phase 
of the complex PQ field, which breaks U(1)pq by its vacuum expectation value 
(VEV). After the end of inflation [2], this phase appears homogenized over 
the universe (supposing that the PQ field is non-zero) with a value 9, which 
is known as the initial misalignment angle. Naturalness suggests that 9 is of 
order unity. This angle remains frozen until the QCD phase transition, where 
the QCD instantons come into play. They break explicitly the PQ symmetry 
to a discrete subgroup ^2 since this symmetry carries QCD anomalies. So, a 
sinusoidal potential for the phase of the PQ field is generated and this phase 
starts oscillating coherently about a minimum of the potential. The resulting 
state resembles pressureless matter consisting of static axions with mass m a ~ 
a qcb/ 'fa, where ^ QCD ~ 200 MeV is the QCD scale. For f a ~ 10 12 GeV, 
the mass of the axion m a ~ 10~ 5 eV. Note that axions, although very light, 
are good candidates for being constituents of the CDM in the universe since 
they are produced at rest. Also, they are very weakly interacting since their 
interactions are suppressed by the axion decay constant f a . 

The relic abundance of axions can be calculated by using the formulae of 
Ref. [TSj, where we take the QCD scale A QC n = 200 MeV and ignore the 
uncertainties for simplicity. We find 



(note that a primordial magnetic helicity, may |14j influence this abundance) . 
So, for natural values of 9 ~ 0.1 and f a ~ 10 12 GeV, axions can contribute 
significantly to CDM, which can even consist solely of axions. 

The main disadvantage of axionic dark matter is that it leads to isocur- 
vature perturbations if the PQ field emerges with non-zero (homogeneous) 
value at the end of inflation. Indeed, during inflation, the angle 9 acquires a 
superhorizon spectrum of perturbations as all the almost massless degrees of 
freedom. At the QCD phase transition, these perturbations turn into isocurva- 
ture perturbations in the axion energy density, which means that the partial 
curvature perturbation in axions is different than the one in photons. The 
recent results of WMAP P put stringent bounds EH El on the possible 
isocurvature perturbation. So, a large axion contribution to CDM is disfa- 
vored in models where the inflationary scale is superheavy (i.e. of the order 
of the SUSY GUT scale) and the PQ field is non-zero at the end of inflation. 

We now wish to turn to the discussion of the main SUSY candidates for 
dark matter: the lightest neutralino x, the axino a and the gravitino G. We 
will consider them mainly within the simplest SUSY framework, which is the 
MSSM. It is, thus, important to first outline the basics of MSSM. 




(1) 
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3 Salient Features of MSSM 

We consider the MSSM embedded in some general SUSY GUT model. We fur- 
ther assume that the GUT gauge group breaking down to the standard model 
(SM) gauge group Gsm occurs in one step at a scale Mqut ~ 10 16 GeV, where 
the gauge coupling constants of strong, weak, and electromagnetic interactions 
unify. Ignoring the Yukawa couplings of the first and second generation, the 
effective superpotential below Mgut is 

W - e t3 (-h t H* 2 qit c + h b H\qib c + h T H\V z r c - ^H{H^), (2) 

where qs = (t, b) and I3 = (v T , r) are the quark and lepton SU(2)l doublet left 
handed supcrfields of the third generation and t c , b c , and r c the corresponding 
SU(2)l singlets. Also, Hi, H 2 are the electroweak Higgs superfields and e the 
2x2 antisymmetric matrix with ei2 = +1. The gravity-mediated soft SUSY- 
breaking terms in the scalar potential are given by 



Koft = ^2 m2 ab ( t > *a ( t>b + 



a.b 



e t] {-A t h t W 2 igt c + A h h b H\qib c + A T h T H\l{f c - BnH\Hi) + h.c.J , (3) 

where the sum is taken over all the complex scalar fields <j> a and tildes denote 
superpartners. The soft gaugino mass terms in the Lagrangian are 

£ g augi„o = \ (m\BB + M 2 J2 W r W r + M 3 ^ ~g a ~g a + h.c. J , (4) 

\ r=l a=l / 

where B, W r and g a are the bino, winos and gluinos respectively. 

The Lagrangian of MSSM is invariant under a discrete Z2 matter parity 
symmetry under which all "matter" (i.e. quark and lepton) superfields change 
sign. Combining this symmetry with the Z2 fermion number symmetry under 
which all fermions change sign, we obtain the discrete Z 2 R-parity symmetry 
under which all SM particles are even, while all sparticles are odd. By virtue of 
R-parity conservation, the LSP is stable and, thus, can contribute to the CDM 
in the universe. It is important to note that matter parity is vital for MSSM 
to avoid baryon- and lepton-number-violating renormalizable couplings in the 
superpotential, which would lead to highly undesirable phenomena such as 
very fast proton decay. So, the possibility of having the LSP as CDM candidate 
is not put in by hand, but arises naturally from the very structure of MSSM. 

The SUSY-breaking parameters m a b, A t ,Ab, A T ,B, and Mi (i = 1, 2, 3) are 
all of the order of the soft SUSY-breaking scale Msusy ~ 1 TeV, but are oth- 
erwise unrelated in the general case. However, if we assume that soft SUSY 
breaking is mediated by minimal supergravity (mSUGRA), i.e. supergrav- 
ity with minimal Kahlcr potential, we obtain soft terms which are universal 
"asymptotically" (i.e. at Mq UT ). In particular, we obtain a common scalar 
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mass Too , a common trilinear scalar coupling A , and a common gaugino mass 
M 1 / 2 - The MSSM supplemented by universal boundary conditions is known as 
constrained MSSM (CMSSM) [HJ. It is true that mSUGRA implies two more 
asymptotic relations: Bq — Aq — mo and too = to 3 / 2 , where Bq — B(Mgvt) 
and TO3/2 is the (asymptotic) gravitino mass. These extra conditions are usu- 
ally not included in the CMSSM. Imposing them, we get the so-called very 
CMSSM JHjj which is a very restrictive version of MSSM and will not be 
considered in these lectures. 

The CMSSM can be further restricted by imposing asymptotic Yukawa 
unification (YU) [201, i.e. the equality of all three Yukawa coupling constants 
of the third family at Mqut: 

h t (M GVT ) = h b (M GVT ) = h T (M GVT ) ee h . (5) 

Exact YU, which makes the CMSSM considerably more predictive, can be 
obtained in GUTs based on a gauge group such as SO(10) or E 6 under which 
all the particles of one family belong to a single representation with the addi- 
tional requirement that the masses of the third family fcrmions arise primarily 
from their unique Yukawa coupling to a single supcrfield representation which 
predominantly contains the electroweak Higgs superfields. It should be noted 
that exact YU in the CMSSM leads to unacceptable values for the bottom- 
quark mass TOf, and, thus, must be corrected in order to become consistent 
with observations. We will ignore this problem for the moment, but we will 
return to it in Sec. [7] 

Now, we assume that our effective theory below M GU t is the CMSSM 
with YU. This theory depends on the following parameters (fj,o = /i(A/gut)): 

m , M1/2, A , B , a G UT, Af G uT, h , tan/3, 

where acuT = ^gut/^ 71 " with <?gut being the GUT gauge coupling constant 
and tan/3 ee (H2) / (Hi) is the ratio of the two electroweak VEVs. The pa- 
rameters o:gut and Mqut are evaluated consistently with the experimental 
values of the electromagnetic and strong fine-structure constants a om and a s , 
and the sine-squared of the Weinberg angle sin 2 9w at Mz- To this end, we 
integrate |21| numerically the renormalization group equations (RGEs) for the 
MSSM at two loops in the gauge and Yukawa coupling constants from Mqut 
down to a common but variable [22] SUSY threshold Msusy = ^/ m t 1 m t 2 (*l,2 
are the stop-quark mass eigenstates). From Msusy to Mz, the RGEs of the 
non-SUSY SM are used. The set of RGEs needed for our computation can be 
found in many references (see e.g. Ref. [2!]). We take a s (Mz) = 0.12 ± 0.001 
which, as it turns out, leads to gauge coupling unification at Mqut with 
an accuracy better than 0.1%. So, we can assume exact unification once the 
appropriate SUSY particle thresholds are taken into account. 

The unified third generation Yukawa coupling constant ho at A/qut 
and the value of tan (3 at Afgusy are estimated using the experimental in- 
puts for the top-quark mass m t (m t ) = 166 GeV and the r-lepton mass 
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m T (Mz) — 1.746 GeV. Our integration procedure of the RGEs relies [22 
on iterative runs of these equations from M GU t to low energies and back for 
every set of values of the input parameters until agreement with the experi- 
mental data is achieved. The SUSY corrections to m T are taken from Ref. [21] 
and incorporated at Msusy. 

Assuming radiative electroweak symmetry breaking, we can express the 
values of the parameters fj, (up to its sign) and P (or, equivalently, the mass 
mA of the CP-odd neutral Higgs boson A) at Msusy in terms of the other 
input parameters by minimizing the tree-level renormalization group (RG) 
improved potential at Msusy- The resulting conditions are 

m^-m^ tan 2 /? 1 2By 2By 

" tan 2 /?-l 2 Mz > ^ ~ m 2 Hl + m% 2 + 2^ ~ m A ' {&) 

where mu x , m# 2 are the soft SUSY-breaking scalar Higgs masses. We can im- 
prove the accuracy of these conditions by including the full one-loop radiative 
corrections to the potential from Ref. [2H a t Msusy- We find that the cor- 
rections to [i and mA from the full one-loop effective potential are minimized 
12211201 by our choice of Msusy- So, a much better accuracy is achieved by 
using this variable SUSY threshold rather than a fixed one. Furthermore, we 
include in our calculation the two-loop radiative corrections to the masses rrih 
and ran of the CP-even neutral Higgs bosons h and P. These corrections are 
particularly important for the mass of the lightest CP-even neutral Higgs bo- 
son h. Finally, the SUSY corrections to mj, are also included at Msusy using 
the relevant formulae of Ref. (21) • As already mentioned, the predicted value 
of the bottom-quark mass is not compatible with experiment. However, we 
will ignore this problem for the moment. The sign of \i is taken to be positive, 
since the fi < case is excluded because it leads [23 [2H1 to a neutralino relic 
abundance which is well above unity, thereby overdosing the universe, for all 
m^'s permitted by b — > sj. We are left with m , Mjy 2 and A as free input 
parameters. 

The LSP is the lightest neutralino \. The mass matrix for the four neu- 
tralinos is 

/ Mi —Mzsw cos M z sw sin (3 \ 

M 2 M z cw cos f3 -M z cw sin (3 



-M z sw cos (3 Mzcw cos f3 —fx 

\ Mzsw sin j3 —Mzcw sin (3 —fx / 



(7) 



in the {—iB, — 1W3, Pi, P2) basis. Here, sw = sinOw, cw = cosOw, and Mi, 
M2 are the mass parameters of P, W3 in Eq. Q). In CMSSM, the lightest 
neutralino turns out to be an almost pure bino P. 

The LSPs are stable due to the presence of the unbroken R-parity, but 
can annihilate in pairs since this symmetry is a discrete Z 2 symmetry. This 
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reduces their relic abundance in the universe. If there exist sparticles with 
masses close to the mass of the LSP, their coannihilation J23J with the LSP 
leads to a further reduction of the LSP relic abundance. It should be noted that 
the number density of these sparticles is not Boltzmann suppressed relative to 
the LSP number density. They eventually decay yielding an equal number of 
LSPs and, thus, contributing to the relic abundance of the LSPs. Of particular 
importance is the next-to-LSP (NLSP), which, in CMSSM, turns out to be 
the lightest stau mass eigenstate T2 . Its mass is obtained by diagonalizing the 
stau mass-squared matrix 

/ m* + + M%(-± + s$y) cos 2/3 m T ( A T - fi tan j3) \ 

(8) 

\ m T (A T — [i tan /3) m 2 . + m? R — Mfs^ cos 2/3 / 

in the gauge basis (tl, tr). Here, m.f L(R] is the soft SUSY-breaking mass of 
the left [right] handed stau tl[RJ an d w r the tau-lepton mass. The stau mass 
eigenstates are 




where Of is the tl — tr mixing angle. 

The large values of b and r Yukawa coupling constants implied by YU cause 
soft SUSY-breaking masses of the third generation squarks and sleptons to run 
(at low energies) to lower physical values than the corresponding masses of the 
first and second generation. Furthermore, the large values of tan/3 implied by 
YU lead to large off-diagonal mixings in the sbottom and stau mass-squared 
matrices. These effects reduce further the physical mass of the lightest stau, 
which is the NLSP. Another effect of the large values of the b and r Yukawa 
coupling constants is the reduction of the mass wia of the CP-odd neutral 
Higgs boson A and, consequently, the other Higgs boson masses to smaller 
values. 



4 Neutralino Relic Abundance 

We now turn to the calculation of the cosmological relic abundance of the 
lightest neutralino x (almost pure B) in the CMSSM with YU according to the 
standard cosmological scenario (for non-standard scenaria, see e.g. Ref. [3D]). 
In general, all sparticles contribute to (l^fi 2 , since they eventually turn into 
LSPs, and all the (co)annihilation processes must be considered. The most 
important contributions, however, come from the LSP and the NLSP. So, in 
the case of the CMSSM, we should concentrate on x (LSP) and f 2 (NLSP) 
and consider the coannihilation of x with fi and f|. The important role of 
the coannihilation of the LSP with sparticles carrying masses close to its mass 
in the calculation of the LSP relic abundance has been pointed out by many 
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authors (see e.g. Refs. [21 [Ml EH EH E3|)- Here, we will use the method of 
Ref. [221 , which was also used in Ref. Note that our analysis can be 

readily applied to any MSSM scheme where the LSP and NLSP are the bino 
and stau respectively. In particular, it applies to the CMSSM without YU, 
where we have tan/3 as an extra free input parameter. 

The relevant quantity, in our case, is the total number density 

n = n% + n? a +n?*, (10) 

since the T2's and f^'s decay into x's after freeze-out. At cosmic temperatures 
relevant for freeze-out, the scattering rates of these (non-relativistic) sparticles 
off particles in the thermal bath are much faster than their annihilation rates 
since the (relativistic) particles in the bath are considerably more abundant. 
Consequently, the number densities Ui (i = x, f%, 7$) are proportional to their 
equilibrium values n° q to a good approximation, i.e. rii/n ~ n^ q /n eq = fj. The 
Boltzmann equation (see e.g. Ref. [31]) is then written as 

^ = -3Hn - (a cS v)(n 2 - (n c *) 2 ), (11) 
at 

where H is the Hubble parameter, v is the "relative velocity" of the annihi- 
lating particles, (• • •) denotes thermal averaging and a e g is the effective cross 
section defined by 

Q"eff = y^^ijUTj (12) 

with (7^ being the total cross section for particle i to annihilate with particle 
j averaged over initial spin states. In our case, a e g takes the following form 

o-cff = a xx r x r x + 4cr xT2 r x r f2 + 2(o> 2f2 + o-fzf* )r f2 r f2 . (13) 
For fj, we use the non-relativistic approximation 

ri{x) = , (14) 

9eS 

g eS (x) = Y / 9i(l + A i )h- A '*, Ai = mi ~ m * . (15) 

• m x 

Here gi = 2, 1, 1 (i = x, f 2 , f| ) is the number of degrees of freedom of the i-th 
particle with mass rrii and x — rn^/T with T being the photon temperature. 

Using Boltzmann equation (which is depicted in Eq. (| 1 If) . we can calculate 
the relic abundance of the LSP at the present cosmic time. It has been found 

El to be given by 

n , 2 1.07 x 10 9 GeV" 1 

%/i 2 fa j (16) 

gi M P Xp X a e {£ 
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with 

poo 

(Toff = xf I (a c gv)x~ 2 dx. (17) 



J XF 

Here Mp ~ 1.22 x 10 19 GeV is the Planck scale, g» ~ 81 is the effective number 
of massless degrees of freedom at freeze-out j3J and xf = m^/Tp with Tp 
being the freeze-out photon temperature calculated by solving iteratively the 
equation [511155] 

0.038 ges(xF) M P (c + 2) cm^ (a e sv)(xF) . . 

x.f = in j — i • l 18 J 

9* X F 

The constant c is chosen to be equal to 1/2 [35] ■ The freeze-out temperatures 
which we obtain here are of the order of m^/25 and, thus, our non-relativistic 
approximation (see Eq. (|14[0 is a posteriori justified. 

Away from s-channel poles and final-state thresholds, the quantities o~iji) 
are well approximated by applying the non-relativistic Taylor expansion up 
to second order in the relative velocity v. 

a tj v = a i0 + b i0 v 2 . (19) 

Actually, this corresponds to an expansion in s and p waves. The thermally 
averaged cross sections are then easily calculated 

x% f°° xv 2 b- 

(aijv)(x) = — -=. \ dvv 2 (a ij v)e~ 2 ^~ = aij + 6— . (20) 
^v 71 " Jo x 

Here, we approximated the masses of the incoming particles by the neutralino 
mass, i.e. rrii = rrij — m^. The reduced mass of the incoming particles is then 
equal to m^/2. We also thermally averaged over the relative velocity rather 
than the separate velocities of the incoming particles, which would be more 
accurate. Using Eqs. lJT7jl. and (J2DJ, one obtains 

^cff = y^X a (ij) a ij + = y^A(ij)> ( 21 ) 

(ij) to') 
where we sum over (ij) = (xx), (x^), and {f 2 fP) with a ( ,> = a f2 f 2 +af 2 f„* , 

b f2f (*) = b f 2 f2 + b f 2 f^ and a (ij)> s iven b y 

r°° dx f°° dx 

a (ij) - c (.ij) x F / —n{x)rj[x), /3 ( y) = 6c {lj) x F / —n(x)rj(x). (22) 

J xf J xf 

Here c (y) = 1, 4, 2 for (ij) = (xx), (x^), and (faf^* 5 ) respectively. 

It should be emphasized that, near s-channel poles or final-state thresholds, 
the Taylor expansion in Eq. (|19fl fails |2H|3Sj badly and, thus, the thermal av- 
erage in Eq. (|20|l has to be calculated accurately by numerical methods. Also, 
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Table 1. Feynman Diagrams 



lllllldl OldLt! 


filial oiaLes 


Diagrams 


XX 


J J 

hh, hH, HH, HA, AA, ZA, ZZ 
hA, hZ, HZ 
H+H~, W+W~ 
W ± H^ 


s(h,H,A,Z), t(J), u{J) 
s(h,H), t( X ), U(x) 
s(A,Z), t(x), u(x) 
s(h,H,Z), t(x±), uCx*) 
s(h,H,A), t( x ±), u( x ±) 


XJi 


rh, tH, tZ 
tA 


s(r), t(ri, 2 ) 

s(t), t(n) 

S(T), t(T 2 ) 


T2T2 


TT 


*(x). <x) 


T2T2 


hh, hH, HH, ZZ 
AA 

H+H~, W+W~ 

77, l z 
it, bb 

TT 

uu, dd, ee 


s(h,H), t(f 1)2 ), «(n,a), c 
s(h,H), i(fi), u(fi), c 
s(h,H,j,Z), t(v T ), c 
t{f 2 ), u(f 2 ), c 
s(h,Hn,Z) 
s(h,H,j,Z), t(x) 



for better accuracy, we should use fully relativistic formulae instead of the 
non-relativistic expressions in Eqs. lfl"3|l . (fl*4"|) . and J20J- Finally, in Eq. (j20J, 

we must take the thermal average over the two initial particle velocities m 
and Vj separately and not just over their relative velocity v. The masses of 
the incoming particles should also be taken different m, ^ nij . After all these 
improvements, Eq. (|20|) takes [37] the form 



1 



x lJ 1 2m 2 i m 2 j TK 2 (^) K 2 (^f) j {mi+mj) 
where K„ are Bessel functions, s the usual Mandelstam variable 



(23) 



s 



m 2 (mf — m 2 ) 2 



^) = i-^^ + ^4^' (24) 



and 



" sW = ») / p$FI^ (2 ' )4i4fe + » - p " )l ^ 12 (25) 

with p', p", E' , E" being the 3-momenta and energies of the outgoing particles 
and \Tij 1 2 the squared transition matrix element summed over final-state spins 
and averaged over initial-state spins. Summation over all final states is implied. 

The relevant final states and Feynman diagrams for % — f 2 (co) annihilation 
are listed in Tabled The exchanged particles are indicated for each pair of 
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initial and final states. The symbols s(x, y, ...), t(x, y, ...), and u(x, y, ...) denote 
tree-level graphs in which the particles x, y, ... are exchanged in the s-, t-, and 
u-channel respectively. The symbol c stands for "contact" diagrams with all 
four external legs meeting at a vertex. The charged Higgs bosons are denoted 
as H^, while / stands for all the matter fermions (quarks and leptons) and e, 
u, and d represent the first and second generation charged leptons, up-, and 
down-type quarks respectively. The bars denote the anti-fermions, are the 
charginos, and v T is the superpartner of the r-neutrino. We have included all 
possible x — X annihilation processes (see e.g. Ref. |HB]), but only the most 
important x~^2, T2, and ti— f| coannihilation processes from Refs. |21II39| 
(for a complete list see e.g. Ref. 0U|)j which are though adequate for giving 
accurate results for all values of tan/3, including the large ones. Some of the 
diagrams listed here have not been considered in previous works 32 , 33 with 
small tan/3. 

The x ~ X annihilation via an A- or _ff-pole exchange in the s-channel can 
be 0J ver y important especially in the CMSSM with large tan/3. As tan/3 
increases, the Higgs boson masses vtia and m# decrease due to the fact that 
hi, increases and, thus, its influence on the RG running of these masses is 
enhanced. As a consequence, ttla and m# approach 2m^ and the neutralino 
pair annihilation via an A- or iJ-pole exchange in the s-channel is resonantly 
enhanced. The contribution from the H pole is p-wave suppressed as one 
can show |31| using CP invariance (recall that the p wave is suppressed by 
xf ~ 25). Therefore, the dominant contribution originates from the A pole 
with the dominant decay mode being the one to bb since, for large tan/3, 
the Abb coupling is enhanced. We find 02] that there exists a region in the 
parameter space of the CMSSM corresponding to large values of tan/3 where 
the x ~ X annihilation via an A pole reduces drastically the relic neutralino 
abundance and, thus, makes it possible to satisfy the WMAP constraint on 
CDM (note that, generically, H^h 2 comes out too large). 

As we already mentioned, near the A pole, the partial wave (or Taylor) 
expansion in Eqs. (|T§)) and (|2*0l fails 12211201 badly. So, the thermal averaging 
must by performed exactly using numerical methods and employing the for- 
mulae in Eqs. (J2HJ), (J2U, and (|2"5|l . In order to achieve good accuracy, it is also 
important to include the one-loop QCD corrections 021 to the decay width of 
the A particle entering its propagator as well as to the quark masses. 

Another phenomenon which helps reducing drastically ft^h? and, thus, 
satisfying the CDM constraint is strong x ~ f i coannihilation |^ 1321 
which operates when nif 2 gets close to m^. This yields [321 G2] a relatively 
narrow allowed region in the m — M1/2 plane (for fixed A and tan/3), which 
stretches just above the excluded region where the LSP is the f 2 . 

There exists 021 also a "bulk" region at m ~ M x /2 ~ few x 100 GeV 
which is allowed by CDM considerations. The (co)annihilation is enhanced 
in this region due to the low values of the various sparticle masses. However, 
this region is, generally, excluded by other phenomcnological constraints (see 
Sec. I7.4JI . So, the A-pole annihilation of neutralinos and the x ~ ^2 coannihi- 
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lation are the two basic available mechanisms for obtaining acceptable values 
for the neutralino relic abundance in the CMSSM. 

There are publicly available codes such as the micrOMEGAs 03] or the 
DarkSUSY 45 for the calculation of fi^h? in MSSM which, among other im- 
provements, include all the relevant (co)annihilation channels between all the 
sparticles (neutralinos, charginos, squarks, sleptons, gluinos), use exact tree- 
level cross sections, calculate accurately and relativistically the thermal aver- 
ages, treat poles and final-state thresholds properly, integrate the Boltzmann 
equation numerically, and include the one-loop QCD corrections to the de- 
cay widths of the Higgs particles and the fermion masses. These codes apply 
to any composition of the neutralino and also include other phenomenolog- 
ical constraints such as the accelerator bounds on certain (s)particle masses 
and the bounds on the anomalous magnetic moment of the muon and the 
branching ration of the process b — > 57 (see Sec. I7.4|) . 

5 Axinos 

Another SUSY particle that could account for the CDM in the universe is 013 
(see also Ref. 07J) the axino a. This particle, which is the superpartner of the 
axion field, is a neutral Majorana chiral fermion with negative R-parity. Its 
mass nia is |48j strongly model-dependent and can be anywhere in the range 
1 eV — MgusY- In the limit of unbroken SUSY, the axino mass is obviously 
equal to the axion mass, which is tiny. Soft SUSY breaking, however, generates 
suppressed corrections to ma via non-renormalizable operators of dimension 
five or higher. So, the corrected mass is at most of order Afg UgY //a ~ 1 keV 
(note that no dimension-four soft mass term is allowed for the axino since 
this particle is a chiral fermion). In specific SUSY models, there also exist 
one- loop contributions to mg, which are typically ^ Mstjsy- When the axion 
is a linear combination of the phases of more than one superfields, we can even 
have tree-level contributions to the axino mass which can easily be as large 
as Msusy- In conclusion, ma is basically a free parameter ranging between 
1 eV and MsuSY' This means that the axino can easily be the LSP in SUSY 
models. 

The axino couplings are suppressed by f a with the most important of them 
being the dimension-five axino (a)-gaugino (A)-gauge boson (A) Lagrangian 
coupling: 

£aXA = ^ a l^ Y °75[7", 71^* + ^frr5[7", TVifc, (26) 

where B and B are, respectively, the gauge boson and gaugino corresponding 
to U(l)y, F b and g b the gluon and gluino fields, ay and a s the U(l)y and 
strong fine-structure constants, and C a yy a model-dependent coefficient of 
order unity. 
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Inflation dilutes utterly any pre-existing axinos, which, after reheating, 
are not in thermal equilibrium with the thermal bath because of their very 
weak couplings (suppressed by /„). They can, however, be thermally produced 
from the bath by 2-body scattering processes or the decay of (s)particles. The 
so-produced axinos are initially relativistic, but out of thermal equilibrium. 
This thermal production (TP) of axinos is @B] predominantly due to 2-body 
scattering processes of strongly interacting particles (because of the relative 
strength of strong interactions) involving the agF coupling in Eq. I|26(l . Such 
processes are 

g + g 5 + g, g + g^a + g, g + q -> a + q, g + q^a + q, 
q + q — > a + g, g + g — > a + g, g + q — > a + q, g + q — > a + q, 

q + q -> a + g, q + q^ a + g, (27) 

where gluons and quarks are denoted by g and q respectively. There exists 46 
also TP of axinos from the decay of thermal gluinos (g — > a + g) or thermal 
neutralinos (x a+"/ \or Z]). The latter proceeds through the dimension- five 
Lagrangian coupling aBB in Eq. (|26|l provided that the neutralino possesses 
an appreciable bino component. These two decay processes are important only 
for reheat temperatures T T of the order of the gluino mass rrig or the neutralino 
mass respectively. 

There is also non-thermal production (NTP) of axinos resulting from the 
decays of sparticles which are out of thermal equilibrium. Indeed, due to 
the suppressed couplings of the axino, the sparticles first decay to the light- 
est ordinary sparticle (LOSP), i.e. the lightest sparticle with non-trivial SM 
quantum numbers, which is the NLSP in this case. The LOSPs then freeze 
out of thermal equilibrium and eventually decay into axinos. 

If the LOSP happens to be the lightest neutralino, the relevant decay pro- 
cess is X ~~ * 5 + 7 [ or Z] through the coupling aBB in Eq. 1)26(1 provided 
that x h as a B component. If, alternatively, the LOSP is the lightest stau mass 
eigenstate, the decay process for the NTP of axinos is T2 — > r + a via the 
one-loop Feynman diagrams in Fig. ^ which contain the effective vertex xfll 
[or xaZ] from the coupling aBB in Eq. (|26l) . In the decay of x, 7's and qq pairs 
are produced. The latter originate from virtual 7 and Z, or real Z exchange 
and lead to hadronic showers. In the T2 case, the resulting t decays immedi- 
ately into light mesons yielding again hadronic showers. The electromagnetic 
and hadronic showers emerging from the LOSP decay in both cases, if they 
are generated after big bang nucleosynthesis (BBN), can cause destruction 
and/or overproduction of some of the light elements, thereby jeopardizing the 
successful predictions of BBN. This implies some constraints on the param- 
eters of the model which, in the present case where the axino is the LSP, 
come basically from the hadronic showers alone due to the relatively short 
LOSP lifetime. In the case of a neutralino LOSP, we obtain 0S| the bound 
rria > 360 MeV for low values of the neutralino mass (< 60 GeV), but no 
bound on the axino mass is obtained for higher values of (> 150 GeV). 
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Fig. 2. The one-loop diagrams for the decay q — > q + a. 



We must further impose the following constraints: (a) the predicted axino 
abundance J7g/i 2 should lie in the 95% c.l. range for the CDM abundance 
in the universe derived by the WMAP satellite [I], (b) both the TP and 
NTP axinos must become non-relativistic before matter domination so as to 
contribute to CDM, and (c) the NTP axinos should not contribute too much 
relativistic energy density during BBN since this can destroy its successful 
predictions. For both \ or T2 LOSP, the requirements (6) and (c) imply that 
ma Z 100 keV or, equivalently, T r ;$ 5x 10 6 GeV. For large values of the reheat 
temperature (T r > 10 4 GeV), TP of axinos is more efficient than NTP and 
the cosmologically favored region in parameter space where the requirement 
(a) holds is quite narrow. For smaller T r 's, NTP dominates yielding a much 
wider favored region with ma > 10 McV. The upper bound on ma increases 
as T r decreases towards m-f. For m§ <C ms, TP of axinos via the process 
q — > q + a becomes |50j very efficient leading to a reduction of the upper limit 
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on T r . As a result, the cosmologically favored region from NTP is reduced in 
this case. The Feynman diagrams for the process q — > q + a are depicted in 
Fig. |21 The restrictions on the ma — T r plane from axino CDM considerations 
are presented in Fig. [3] 




-^qO I I I I I I 

icr 4 10" 2 io° io 2 

m~ a (GeV) 

Fig. 3. The restrictions on the ma — T r plane from axino CDM considerations 
for x = B, m s = 100 GeV, m- g = m 9 = 1 TeV, and f a = 10 11 GeV. The solid 
almost diagonal line corresponds to f2j F h 2 w 0.13, where Hj p h 2 is the TP axino 
abundance. So, the area above this line is cosmologically excluded. The narrow 
shaded area just below the thin part of this line for ma < 10 MeV is cosmologically 
favored by TP. The hatched areas are favored by NTP. For m, <C m§, the solid line 
is replaced by the dashed one, whose position is strongly dependent on the actual 
values of m ? -, m§ and is only indicative here. The area favored by NTP is then 
limited only to the "back-hatched" region which lies below the dashed line. 

We find 0H| that, for the CMSSM, with appropriate choices of ma and 
T r , almost any pair of values for mo and M 1 / 2 can be allowed. This holds for 
both x or T2 as LOSP. However, the required T r 's for achieving the WMAP 
bound on CDM turn out to be quite low few x 100 GeV). 
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6 Gravitinos 

It has been proposed that CDM could also consist of gravitinos. The 

gravitino G is the superpartner of the graviton and has negative R-parity. It 
can be the LSP in many cases and, thus, contribute to CDM. In the very 
CMSSM, its mass rriQ is fixed by the asymptotic condition m 3 / 2 = tuq. In the 
general CMSSM, however, it is a free parameter ranging between 100 GeV 
and 1 TeV. It can, thus, very easily be the LSP in this case. 

The couplings of the gravitino are suppressed by the Planck scale. The 
most important of them are given by the dimension-five Lagrangian terms 

c = _J_x>„^ tW - -J—x^VVyV/i 

V Imp v 2mp 

-^^b u ,i p ]r~^ p , (28) 

where ip^ denotes the gravitino field, (ft are the complex scalar fields, ip l 
are the corresponding chiral fermion fields, A a are the gaugino fields, mp ~ 
2.44 x 10 18 GeV is the reduced Planck scale, and T> v denotes the covariant 
derivative. From these Lagrangian terms, we obtain scalar-fermion-gravitino 
vertices (0/G) such as qqG, UG, and HHG, as well as gaugino-gauge boson- 
gravitino vertices (XFG) such as ggG and BBG (in this section, I and H 
represent any lepton and Higgs boson respectively). 

The gravitinos are thermally produced after reheating by 2 — > 2 scattering 
processes involving the above vertices. Such processes are 152) 

G + q, 
G + q, 

(29) 

There is jS3 ESj also NTP of gravitinos via the decay of the NLSP. For 
neutralino NLSP, the relevant decay processes are x ~~ * G + 7 [or Z] from the 
XFG coupling and x ~~ * G + H from the HHG coupling. In the case of f% 
NLSP, the relevant decay process is T2 — > t+G from the vertex UG. There is an 
important difference between the NTP of gravitinos and axinos. In the former 
case, the NLSP has a large lifetime (up to about 10 8 sec). Consequently, it 
gives rise mostly to electromagnetic, but also to hadronic showers well after 
BBN. The electromagnetic showers cause destruction of some light elements 
(D, 4 He, 7 Li) and/or overproduction of 3 He and 6 Li, thereby disturbing BBN. 
The hadronic showers can also disturb BBN. The overall resulting constraint 
is |54| very strong allowing only limited regions of the parameter space of the 
CMSSM lying exclusively in the range where the NLSP is the f^. Moreover, 
in these allowed regions, the NTP of gravitinos is not efficient enough to 
account for the observed CDM abundance for M 1 / 2 ~ 6 TeV. However, we 
can compensate for the inefficiency of NTP by raising T r to enhance the TP 
of G's. The relic gravitino abundance from TP, for ttlq <C to§, is 
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where Wg(/i) is the running gluino mass (for the general formula, see Ref. |56p. 
7 Yukawa Quasi-Unification 

As already said in Sec. |31 exact YU in the framework of the CMSSM leads to 
wrong values for nib and, thus, must be corrected. We will now present a model 
which naturally solves |25] (see also Refs. |57U58 |) this nib problem and discuss 
the restrictions on its parameter space implied by CDM considerations and 
other phenomenological constraints. Exact YU can be achieved by embedding 
the MSSM in a SUSY GUT model with a gauge group containing SU(4) c and 
SU(2) R . Indeed, assuming that the electroweak Higgs superfields Hi, Hi and 
the third family right handed quark superfields t c , b c form SU(2) R doublets, 
we obtain |59| the asymptotic Yukawa coupling relation h t = hb and, hence, 
large tan/3 ~ m t /mb. Moreover, if the third generation quark and lepton 
SU(2) L doublets [singlets] q 3 and l 3 [b c and r c ] form a SU(4) c 4-plet [4-plet] 
and the Higgs doublet Hi which couples to them is a SU(4) c singlet, we get 
hb = h T and the asymptotic relation nib — m, T follows. The simplest GUT 
gauge group which contains both SU(4) C and SU(2) R is the Pati-Salam (PS) 
group Gps = SU(4) C x SU(2) L x SU(2) R and we will use it here. 

As mentioned, applying YU in the context of the CMSSM and given the 
experimental values of the top-quark and tau-lepton masses (which naturally 
restrict tan/3 ~ 50), the resulting value of the &-quark mass turns out to be 
unacceptable. This is due to the fact that, in the large tan/3 regime, the tree- 
level &-quark mass receives sizeable SUSY corrections |2H EHH El ES2I (about 
20%), which have the sign of \i (with the standard sign convention [HSj) and 
drive, for u, > [<] 0, the corrected &-quark mass at Mz, rrib(Mz), well above 
[somewhat below] its 95% c.l. experimental range 

2.684 GeV < m b (M z ) < 3.092 GeV with a s {M z ) = 0.1185. (31) 

This is derived by appropriately |39| evolving the corresponding range of 
mb(rrib) in the MS scheme (i.e. 3.95 — 4.55 GeV) up to Mz in accordance 
with Ref. |64j . We see that, for both signs of fi, YU leads to an unacceptable 
6-quark mass with the (i < case being less disfavored. 

A way out of this mb problem can be found [3^1 (see also Refs. [571 155] ]) 
without having to abandon the CMSSM (in contrast to the usual strategy 
|ri2l 1651 lr>6l Ifi7| ) or YU altogether. We can rather modestly correct YU by 
including an extra SU(4) C non-singlet Higgs superfield with Yukawa couplings 
to the quarks and leptons. The Higgs SU(2) L doublets contained in this super- 
field can naturally develop |68| subdominant VEVs and mix with the main 
electroweak doublets, which are assumed to be SU(4) C singlets and form a 
SU(2) R doublet. This mixing can, in general, violate the SU(2) R symmetry. 
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Consequently, the resulting electroweak Higgs doublets H 1 , H 2 do not form a 
SU(2) R doublet and also break the SU(4) C symmetry. The required deviation 
from YU is expected to be more pronounced for /i > 0. Despite this, we will 
study here this case, since the \x < case has been excluded j^H] by combining 
the WMAP restrictions £Q on the CDM in the universe with the experimen- 
tal results [23 on the inclusive branching ratio BR(& — > sj). The same SUSY 
GUT model which, for /i > and universal boundary conditions, remedies 
the nib problem leads to a new version [7T] of shifted hybrid inflation [72] . 
which, as the older version |72) . avoids monopole overproduction at the end 
of inflation, but, in contrast to that version, is based only on renormalizablc 
interactions. 

In Sec. 17.11 we review the construction of a SUSY GUT model which 
naturally and modestly violates YU, yielding an appropriate Yukawa quasi- 
unification condition (YQUC), which is derived in Sec. 17.21 We then outline 
the resulting CMSSM in Sec. 17.31 and introduce the various cosmological and 
phenomenological requirements which restrict its parameter space in Sec. 17.41 
In Sec. 17.51 we delineate the allowed range of parameters. Finally, in Sec. 17.61 
we briefly comment on the new version of shifted hybrid inflation which is 
realized in this model. 

7.1 The PS SUSY GUT Model 

We will take the SUSY GUT model of shifted hybrid inflation [7J] ( see also 
Ref. [SI) as our starting point. It is based on Gps, which is the simplest GUT 
gauge group that can lead to exact YU. The representations under Gps and 
the global charges of the various matter and Higgs superfields contained in 
this model are presented in Table [21 which also contains the extra Higgs su- 
perfields required for accommodating an adequate violation of YU for /x > 
(see below). The matter superfields are Fi and Ff (i = 1, 2, 3), while the elec- 
troweak Higgs doublets belong to the superfield h. So, all the requirements 
for exact YU are fulfilled. The spontaneous breaking of Gps down to Gsm is 
achieved by the superheavy VEVs (~ A/gut) of the right handed neutrino- 
type components of a conjugate pair of Higgs superfields H c , H c . The model 
also contains a gauge singlet S which triggers the breaking of Gps, a SU(4) C 
6-plet G which gives [7J] masses to the right handed down-quark-type com- 
ponents of H c , H c , and a pair of gauge singlets N, N for solving [E] the (x 
problem of the MSSM via a PQ symmetry (for an alternative solution of the 
(U problem, see Ref. [75]). In addition to Gps, the model possesses two global 
U(l) symmetries, namely a R and a PQ symmetry, as well as the discrete mat- 
ter parity symmetry Z 1 ^ ' . Note that global continuous symmetries such as our 
PQ and R symmetry can effectively arise [77] from the rich discrete symmetry 
groups encountered in many compactified string theories (see e.g. Ref. |78|V 
Note that, although the model contains baryon- and lepton-number-violating 
superpotential terms, the proton is E2] practically stable. The baryon 
asymmetry of the universe is generated via the non-thermal realization |7H| of 
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the leptogenesis scenario jSU] (for recent papers on non-thermal leptogenesis, 
see e.g. Ref. EJ). 



Table 2. Superfield Content of the Model 



Superfields 


Representations 


Global 




under Gps 


Charges 






R 


PQ Z™ p 


Matter Superfields 


Fi 


(4,2,1) 


1/2 


-1 1 




(4,1,2) 


1/2 


-1 


Higgs Superfields 


h 


(1,2,2) 





1 


H c 


(4,1,2) 








H c 


(4,1,2) 








S 


(1,1,1) 


1 





G 


(6,1,1) 


1 





N 


(1,1,1) 


1/2 


-1 


N 


(1,1,1) 





1 


Extra Higgs Superfield 


3 




ti 


(15,2,2) 





1 


h' 


(15,2,2) 


1 


-1 




(15,1,3) 








4> 


(15,1,3) 


1 






A moderate violation of exact YU can be naturally accommodated in this 
model by adding a new Higgs superfield h' with Yukawa couplings FF c h'. 
Actually, (15,2,2) is the only representation of Gps, besides (1,2,2), which 
possesses such couplings to the matter superfields. In order to give superheavy 
masses to the color non-singlet components of h', we need to include one more 
Higgs superfield hi with the superpotential coupling h'h' , whose coefficient is 
of the order of Mgut- 

After the breaking of Gps to Gsm, the two color singlet SU(2) L doublets 
h[, h' 2 contained in hi can mix with the corresponding doublets hi, h% in h. 
This is mainly due to the terms h'h' and H c H c h'h. Actually, since 

H C H C = (4, 1, 2) (4, 1, 2) = (15, 1, 1 + 3) + • ■ ■ , 

h'h = (15, 2, 2)(1, 2, 2) = (15, 1, 1 + 3) + • • • , (32) 

there are two independent couplings of the type H c H c h'h (both suppressed 
by the string scale Ms ~ 5 x 10 17 GeV, as they are non-renormalizable) . 



20 George Lazarides 



One of these couplings is between the SU(2) R singlets in H C H C and h'h and 
the other between the SU(2) R triplets in these combinations. So, we obtain 
two bilinear terms h\h\ and h' 2 h 2 with different coefficients, which are sup- 
pressed by Mqut/Ms- These terms together with the terms h\h\ and h' 2 h' 2 
from h'h', which have equal coefficients, generate different mixings between 
hi, h\ and hi, h' 2 . Consequently, the resulting electroweak doublets H\, H 2 
contain SU(4) C violating components suppressed by Mgut /Ms and fail to 
form a SU(2) R doublet by an equally suppressed amount. So, YU is naturally 
and moderately violated. Unfortunately, as it turns out, this violation is not 
adequately large for correcting the bottom-quark mass within the framework 
of the CMSSM with ft > 0. 

In order to allow for a more sizable violation of YU, we further extend the 
model by including the superfield tf> with the coupling <ph'h. To give super- 
heavy masses to the color non-singlets in <p, we introduce one more superfield 
4> with the coupling <jxf), whose coefficient is of order Mgut- 

The superpotential terms (jxp and <pH c H c imply that, after the breaking 
of Gps to Gsm, 4> acquires a VEV of order Mgut- The coupling (fh'h then 
generates SU(2) R violating unsuppressed bilinear terms between the doublets 
in h' and h. These terms can overshadow the corresponding ones from the 
non-renormalizable term H c H c h'h. The resulting SU(2) R violating mixing of 
the doublets in h and b! is then unsuppressed and we can obtain stronger 
violation of YU. 



7.2 The YQUC 

To further analyze the mixing of the doublets in h and h', observe that the 
part of the superpotential corresponding to the symbolic couplings h'h', tfih'h 
is properly written as 

mtr (h'eh' T e) + ptr (h'e(j)h T e) , (33) 

where m is a mass parameter of order Mqut, P is a dimensionless parameter 
of order unity, tr denotes trace taken with respect to the SU(4) c and SU(2) L 
indices, and the superscript T denotes the transpose of a matrix. 

After the breaking of Gps to Gsm, acquires a VEV (</>) ~ Mgut- Sub- 
stituting it by this VEV in the above couplings, we obtain 

tr{h'eh' T e) = \h' 2 + h[ T eh' 2 + ■■■, (34) 

tvih'e^e) = -^ltr(/iW T e) = ^{h'^ eh 2 - hjeh' 2 ), (35) 
V 2 v 2 

where the ellipsis in Eq. I|34|l contains the colored components of h', h' and 
0-3 = diag(l, —1). Inserting Eqs. (|34|) and l|35|) into Eq. (|33|l . we obtain 

mft/ 1 T e(h' 2 — ah 2 ) + m(/i' 1 T + ahj)eh' 2 with a = — ^ ~ . (36) 

V2m 
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So, we get two pairs of superheavy doublets with mass m. They are predom- 
inantly given by 



h'i, 



h' 2 — ah 2 



and 



h'i + ah 1 -, 



ho. 



(37) 



The orthogonal combinations of hi, h\ and h 2 , h' 2 constitute the electroweak 
doublets 



hi - a*h[ h 2 +a*h' 2 
Hi = — and H 2 — 



(38) 



The superheavy doublets in Eq. Q37J1 must have vanishing VEVs, which readily 
implies that (h[) = —a (hi) and (h' 2 ) — a(h 2 )- Equation 138JI then gives 
(Hi) = (1 + |a| 2 ) 1/2 (M, (H?) = (1 + \a\ 2 ) 1/2 (h2l- From the third generation 
Yukawa couplings y^F^hF^ , iy'^F^h' , we obtain 



m t = \yw(h 2 ) +y'^(h' 2 ) \ = 
1 — pa I \/3 



1 + pa/V3 



V33(H 2 



m b = 



-y^(Hi 



1 + \/3pa 



2/33 (Hi 



(39) 
(40) 



where p = 2/33/1/33. From Eqs. l|3l3|) and I^UJ, we see that YU is now replaced 
by the YQUC 

h t : h b : h T = (1 + c) : (1 - c) : (1 + 3c) with < c = pa/Vij < 1. (41) 

For simplicity, we restricted ourselves here to real values of c only which 
lie between zero and unity, although c is, in general, an arbitrary complex 
quantity with |c| ~ 1. 



7.3 The Resulting CMSSM 

Below the GUT scale Mgut, the particle content of our model reduces to this 
of MSSM (modulo SM singlets). We assume universal soft SUSY breaking 
scalar masses mo, gaugino masses M 1( / 2 , and trilinear scalar couplings Aq at 
A/gut- Therefore, the resulting MSSM is the so-called CMSSM [T5] with fx > 
supplemented by the YQUC in Eq. I|41|) . With these initial conditions, we run 
the MSSM RGEs [2] between M GUT and a common variable SUSY threshold 
AfsuSY ( see S ec - 13 determined in consistency with the SUSY spectrum of 
the model. At Msusy, we impose radiative electroweak symmetry breaking, 
evaluate the SUSY spectrum and incorporate the SUSY corrections P4llBTllrj2] 
to the 6-quark and r-lepton masses. Note that the corrections to the r-lepton 
mass (almost 4%) lead [HHl to a small reduction of tan/3. From Msusy to Mz, 
the running of gauge and Yukawa coupling constants is continued using the 
SM RGEs. 
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For presentation purposes, M-^i 2 and m Q can be replaced |21j by the LSP 
mass TOlsp and the relative mass splitting between this particle and the light- 
est stau Af 2 = (rrif 2 — tolsp)/ to lsp (recall that f 2 is the NLSP in this case). 
For simplicity, we restrict this presentation to the Aq = case (for Aq ^ 
see Refs. [3§] EH)- So, our input parameters are TOlsp, A? 2 , c, and tan/3. 

For any given rrib(Mz) in the range in Eq. i|31[) and with fixed m t (m t ) = 
166 GeV and m T (Mz) — 1.746 GeV, we can determine the parameters c and 
tan/3 at Msusy so that the YQUC in Eq. I|41|) is satisfied. We are, thus, left 
with tolsp arid Af 2 as free parameters. 

7.4 Cosmological and Phenomenological Constraints 

Restrictions on the parameters of our model can be derived by imposing a 
number of cosmological and phenomenological requirements (for similar recent 
analyses, see Refs. [SHI El IBS])- These constraints result from 

• CDM Considerations. As discussed in Sec. |3] in the context of the 
CMSSM, the LSP can be the lightest neutralino which is an almost pure bino. 
It naturally arises as a CDM candidate. We require its relic abundance, 
^hsph 2 , not to exceed the 95% c.l. upper bound on the CDM abundance 
derived [Tj by WMAP: 

ficDMft 2 £ 0.13. (42) 

We calculate Sl^pti 2 using micrOMEGAs |44|. which is certainly one of the 
most complete publicly available codes. Among other things, it includes all 
possible coannihilation processes (33J and one-loop QCD corrections 03] to 
the Higgs decay widths and couplings to fermions. 

• Branching Ratio ofb — > sj. Taking into account the experimental results 
of Ref. [S3] on this ratio, BR(b — > sj), and combining (3j| appropriately the 
experimental and theoretical errors involved, we obtain the 95% c.l. range 

1.9 x 1(T 4 < BR(b -► sj) < 4.6 x 1CT 4 . (43) 

Although there exist more recent experimental data [85] on the branching 
ratio of b — > sj, we do not use them here. The reason is that these data do 
not separate the theoretical errors from the experimental ones and, thus, the 
derivation of the 95% c.l. range is quite ambiguous. In any case, the 95% c.l. 
limits obtained in Ref. (HE] on the basis of these latest measurements are not 
terribly different from the ones quoted in Eq. I|43|l . In view of this and the 
fact that, in our case, the restrictions from BR(6 — > sj) are overshadowed 
by other constraints (see Sec. 17. 5f) . we limit ourselves to the older data. We 
compute BR(6 — > sj) by using an updated version of the relevant calculation 
contained in the micrOMEGAs package 03]. In this code, the SM contribution 
is calculated following Ref. jEZ]- The charged Higgs (-H^) contribution is eval- 
uated by including the next-to-leading order (NLO) QCD corrections ]HH] and 
tan /3 enhanced contributions |HH] ■ The dominant SUSY contribution includes 
resummed NLO SUSY QCD corrections [HE] , which hold for large tan /3. 
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• Muon Anomalous Magnetic Moment. The deviation, 8a„, of the mea- 
sured value of a M from its predicted value in the SM, a^ M , can be attributed 
to SUSY contributions, which are calculated by using the micrOMEGAs routine 
|59") . The calculation of a^ M is not yet stabilized mainly because of the insta- 
bility of the hadronic vacuum polarization contribution. According to recent 
calculations (see e.g. Refs. jHOll^), there is still a considerable discrepancy 
between the findings based on the e + e~ annihilation data and the ones based 
on the r-decay data. Taking into account the results of Ref. [00] and the ex- 
perimental measurement of a M reported in Ref. |02], we get the following 95% 
c.l. ranges: 

- 0.53 x 1CT 10 < Saf, < 44.7 x 1(T 10 , e + e~-based; (44) 
-13.6 x 1CT 10 < 8a^ < 28.4 x 1(T 10 , r-based. (45) 

Following the common practice |83|. we adopt the restrictions to parameters 
induced by Eq. (|4~4fl , since Eq. (|43|l is considered as quite oracular, due to poor 
r-decay data. It is true that there exist more recent experimental data [OH] 
on <z M than the ones we considered and more updated estimates of 5a^ than 
the one in Ref. [00] (see e.g. Ref. [05 )• However, only the 95% c.l. upper limit 
on (5a M enters into our analysis here and its new values are not very different 
from the one in Eq. (|44|l . 

• Collider Bounds. Here, as it turns out, the only relevant collider bound 
is the 95% c.l. LEP lower bound [01] on the mass of the lightest CP-even 
neutral Higgs boson h: 

m h > 114.4 GeV. (46) 

The SUSY corrections to the lightest CP-even Higgs boson mass are 
calculated at two loops by using the FeynHiggsFast program |95| included in 
the micrOMEGAs code |4*4*| . 

7.5 The Allowed Parameter Space 

We will now try to delineate the parameter space of our model with /i > 
which is consistent with the constraints in Sec. 17.41 The restrictions on the 
wlsp — Af% plane, for Ao = and the central values of a s (Mz) = 0.1185 and 
mb(Mz) = 2.888 GeV, are indicated in Fig.0]by solid lines, while the upper 
bound on m LS p from Eq. g2J, for m b (M z ) = 2.684 [3.092] GeV, is depicted 
by a dashed [dotted] line. We observe the following: 

• The lower bounds on ttilsp are not so sensitive to the variations of 
m b (M z ). 

• The lower bound on rriLSP from Eq. Q46J) overshadows all the other lower 
bounds on this mass. 

• The upper bound on tolsp from Eq. 142|) is very sensitive to the variations 
of mb(Mz). In particular, one notices the extreme sensitivity of the almost 
vertical part of the corresponding line, where the LSP annihilation via an 
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Fig. 4. The various restrictions on the 771lsp — Af 2 plane for fi > 0, Ao = 0, and 
a s (Mz) = 0.1185. From left to right, the solid lines depict the lower bounds on 
m L sp from 5a„ < 44.7 x KT 10 , BR(6 -> s-y) > 1.9 x 10 -4 , and m h > 114.4 GeV 
and the upper bound on ttilsp from ^lsp^ 2 < 0.13 for mt(Mz) = 2.888 GeV. 
The dashed [dotted] line depicts the upper bound on mLSP from Q^sph 2 < 0.13 for 
m b (M z ) = 2.684 [3.092] GeV. The allowed area for m b (M z ) = 2.888 GeV is shaded. 

A-boson exchange in the s-channel is |96j by far the dominant process, 
since tua, which is smaller than 2?tilsp, is always very close to it as seen 
from Fig. [S] This sensitivity can be understood from Fig. [fJl where tua 
is depicted versus mLSP for various mb(Mz)'s. We see that, as m b (Mz) 
decreases, tua increases and approaches 2?tilsp- The A-pole annihilation 
is then enhanced and X2lsp^ 2 is drastically reduced causing an increase of 
the upper bound on ttilsp- 
• For low Af 2 's, bino-stau coannihilations take over leading to a very 
pronounced reduction of the LSP relic abundance i?LSP^ 2 , thereby en- 
hancing the upper limit on mLSP- So, we obtain the almost horizontal tail 
of the allowed region in Fig. 0] 

For (i > 0, A = 0, a s (M z ) = 0.1185 and m b (M z ) = 2.888 GeV, we find 
the following allowed ranges of parameters: 

176 GeV < m LSP < 615 GeV, < A? 2 < 1.8, 

58 < tan/3 < 59, 0.14 < c < 0.17. (47) 

The splitting between the bottom (or tau) and top Yukawa coupling constants 
Sh= —(hi, — h t )/h t — (h T — h t )/h t = 2c/(l + c) ranges between 0.25 and 0.29. 
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Fig. 5. The mass parameters uia and -Msusy versus to-lsp for various values of Af 2 , 
which are indicated on the curves. We take fi > 0, Ao = 0, rrib(Mz) = 2.888 GeV, 
and a s (M z ) = 0.1185. 




Fig. 6. The mass parameters tua and Msusy as functions of mLSP for H > 0, 
A = 0, A f2 = 1, a 3 (M z ) = 0.1185, and with m b (M z ) = 2.684 GeV (dashed lines), 
3.092 GeV (dotted lines), or 2.888 GeV (solid lines). 
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7.6 The New Shifted Hybrid Inflation 

It is interesting to note that our SUSY GUT model gives rise [7J naturally to 
a modified version of shifted hybrid inflation [72] . Hybrid inflation [57j , which 
is certainly one of the most promising inflationary scenarios, uses two real 
scalars: one which provides the vacuum energy density for driving inflation 
and a second which is the slowly varying field during inflation. This scheme, 
which is naturally incorporated [HE] m SUSY GUTs (for an updated review, 
see Ref. ^99 ), in its standard realization has the following property |f00j : 
if the GUT gauge symmetry breaking predicts topological defects such as 
magnetic monopoles jlUlj . cosmic strings |1(J2) . or domain walls |ll)3j . these 
defects are copiously produced at the end of inflation. In the case of monopoles 
or walls, this leads to a cosmological catastrophe |1U4| . The breaking of the 
Gps symmetry predicts the existence of doubly charged monopoles |105| . So, 
any PS SUSY GUT model incorporating the standard realization of SUSY 
hybrid inflation would be unacceptable. One way to remedy this is to invoke 
thermal inflation jl()7| to dilute the primordial monopoles well after 
their production. Alternatively, we can construct variants of the standard 
SUSY hybrid inflationary scenario such as smooth |100j or shifted |72] hybrid 
inflation which do not suffer from the monopole overproduction problem. In 
the latter scenario, we generate [7_2j a shifted inflationary trajectory so that 
Gps is already broken during inflation. This could be achieved [22] in our 
SUSY GUT model even before the introduction of the extra Higgs superfields, 
but only by utilizing non-renormalizable terms. The inclusion of h! and h' does 
not change this situation. The inclusion of <f> and </>, however, very naturally 
gives rise |71j to a shifted path, but now with renormalizable interactions 
alone. 



8 Conclusions 

We showed that particle physics provides us with a number of candidate par- 
ticles out of which the CDM of the universe can be made. These particles 
are not invented solely for explaining the CDM, but they are naturally there 
in various particle physics models. We discussed in some detail the major 
candidates which are the axion, the lightest neutralino, the axino, and the 
gravitino. The last three particles exist only in SUSY theories and can be 
stable provided that they are the LSP. 

The axion is a pseudo Nambu-Goldstone boson associated with the sponta- 
neous breaking of a PQ symmetry. This is a global anomalous U(l) symmetry 
invoked to solve the strong CP problem. It is, actually, the most natural so- 
lution to this problem which is available at present. The axions are extremely 
light particles and are generated at the QCD phase transition carrying zero 
momentum. We argued that these particles can easily provide the CDM in 
the universe. However, if the PQ field emerges with non-zero value at the end 
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of inflation, they lead to isocurvature perturbations, which, for superheavy 
inflationary scales, are too strong to be compatible with the recent results of 
the WMAP satellite on the CMBR anisotropies. 

The most popular CDM candidate is, certainly, the lightest neutralino 
which is present in all SUSY models and can be the LSP for a wide range of 
parameters. We considered it within the simplest SUSY framework which is 
the MSSM whose salient properties were summarized. We used exclusively the 
constrained version of MSSM which is known as CMSSM and is based on uni- 
versal boundary conditions. In this case, the lightest neutralino is an almost 
pure bino, whereas the NLSP is the lightest stau. We sketched the calculation 
of the neutralino relic abundance in the universe paying particular attention 
not only to the neutralino pair annihilations, but to the neutralino-stau coan- 
nihilations too. It is very important for the accuracy of the calculation to 
treat poles and final-state thresholds properly and include the one-loop QCD 
corrections to the Higgs boson decay widths and the fermion masses. We find 
that two effects help us reduce the neutralino relic abundance and satisfy the 
WMAP constraint on CDM: the resonantly enhanced neutralino pair anni- 
hilation via an A-pole exchange in the s-channel, which appears in the large 
tan (3 regime, and the strong neutralino-stau coannihilation, which is achieved 
when these particles are almost degenerate in mass. 

The axino, which is the SUSY partner of the axion, can also be the LSP 
in many cases since its mass is a strongly mo del- dependent parameter in the 
CMSSM. It is produced thermally by 2-body scattering or decay processes 
in the thermal bath, or non-thermally by the decay of sparticles which are 
already frozen out of thermal equilibrium. For small axino masses, TP is more 
important yielding a very narrow favored region in the parameter space. For 
larger axino masses, however, NTP is more efficient and the favored region 
in the parameter space becomes considerably wider. One finds that, in the 
case of the CMSSM, almost any point on the mo — Mi/ 2 plane can be allowed 
by axino CDM considerations. The required reheat temperatures though are 
quite small (< few x 100 GeV). 

The mass of the gravitino is a practically free parameter in the CMSSM. 
So, the gravitino can easily be the LSP and, in principle, contribute to the 
CDM of the universe. It is produced thermally by 2-body scattering processes 
in the thermal bath as well as non-thermally by the decay of the NLSP, which 
can be either the neutralino or the stau. In contrast to the axino case, however, 
the NLSP can now have quite a long lifetime. The electromagnetic showers 
resulting from the NLSP decay can destroy the successful predictions of BBN. 
So, we obtain strong constraints which allow only very limited regions of the 
parameter space of the CMSSM. As it turns out, NTP in these regions is not 
efficient enough to account for CDM. We can, however, make these regions 
cosmologically favored by raising T r to enhance TP of gravitinos. 

We studied the CMSSM with fi > and A = applying a YQUC which 
originates from a PS SUSY GUT model. This condition yields an adequate 
deviation from YU which allows an acceptable mb(Mz). We, also, imposed 
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the constraints from the CDM in the universe, b — > s"f, 8a „ and m^. We found 
that there exists a wide and natural range of CMSSM parameters which is 
consistent with all the above constraints. The parameter tan j3 ranges between 
about 58 and 59 and the asymptotic splitting between the bottom (or tau) and 
the top Yukawa coupling constants varies in the range 25 — 29% for central val- 
ues of vtibiMz) and a s (Mz). The predicted LSP mass can be as low as about 
176 GeV. Moreover, the model resolves the [i problem of MSSM, predicts 
stable proton, generates the baryon asymmetry of the universe via primordial 
leptogenesis, and gives rise to a new version of shifted hybrid inflation which 
is based solely on renormalizable interactions. 
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